the topological implications of inhomogeneity 



Boudewijn F. Roukema 1 and Vincent Blanlceil 2 and Jan J. 
Ostrowski 1 

1 Toruri Centre for Astronomy, Nicolaus Copernicus University, 
ul. Gagarina 11, 87-100 Toruri, Poland 

2 Departement de Mathematiques, Universite de Strasbourg, 7 
rue Rene Descartes, 67084 Strasbourg cedex, France 

Abstract. The approximate homogeneity of spatial sections of the Universe is well 
supported observationally, but the inhomogeneity of spatial sections is even better 
supported. Here, we consider the implications of inhomogeneity in dust models for 
the connectedness of spatial sections at early times. We consider a Lemaitre-Tolman- 
Bondi (LTB) model designed to match observations, a more general, heuristic model 
motivated by the former, and two specific, global LTB models. We propose that 
the generic class of solutions of the Einstein equations projected back in time from 
the spatial section at the present epoch includes subclasses in which the spatial section 
evolves (with increasing time) smoothly (i) from being disconnected to being connected, 
or (ii) from being simply connected to being multiply connected. We show that (i) 
and (ii) each contain at least one exact solution. These subclasses exist because the 
Einstein equations allow non-simultaneous big bang times. The two types of topological 
evolution occur at post-quantum epochs if the bang time varies by much more than a 
Planck time. Both require physics beyond the Einstein equations. A phenomenological 
outline for modelling the evolution from a simply connected spatial section to a multiply 
connected section is proposed. 
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1. Introduction 

The approximate homogeneity of spatial sections (hypersurfaces) of the Universe 
is well supported observationally. Both the assumption of homogeneity and fact 
of mhomogeneity play an important role in relativistic cosmological models. The 
Friedmann-Lemaitre-Robertson- Walker (FLRW) models 

be, sag 

are solutions 

of the Einstein equations in which the density is constant in any comoving spatial 
section. With the Concordance Model parameters of the metric the FLRW models 
provide reasonably good fits to observational data (faint galaxy number counts [e.g. 
3. jl, gravitational lensing [e.g. 9|, supernovae type la magnitude-redshift relations [e.g. 
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lf|). However, the Universe is clearly inhomogeneous, and the forcing of an FLRW 
model onto the data requires a non-zero "dark energy" parameter Q\, suggesting that 
the latter is most simply interpreted as an artefact of forcing an oversimplified model 
onto the data [e.g. 
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The near-homogeneity is also a key element of the "Horizon Problem" for non- 
inflationary FLRW models: how was it possible for causally disconnected (but spatially 
connected) regions of the spatial section of the Universe to homogenise? In the context 
of dust models with comoving spatial sections, this question implicitly assumes that 
universe models with initially inhomogeneous spatial sections are relativistically valid 
and only have a problem with causal disconnectedness, not with comoving spatial 
disconnectedness. Is this assumption correct? 

Generic inhomogeneous cosmological solutions of the Einstein equations are 
not known. However, the Lemaitre-Tolman-Bondi (LTB) family of exact solutions 



15l . 1 161 . Il7j is known. These solutions consist of exact solutions to the Einstein equations 
that are radially inhomogeneous and spherically symmetric with respect to an origin. In 
analogy with the way that the FLRW model is interpreted to apply to a 3-dimensionally 
averaged spatial section, an LTB solution can be interpreted to apply to a spatial 
section that has been averaged over every infinitesimally thin spherical shell, i.e. 2- 
dimensionally, prior to solving the Einstein equations. However, the Einstein equations 
do not imply their averaged equivalent: G(g) = 87rT(p) G((g)) = 87rT((p)j|| [e.g. 



181 ]. That is, although homogeneity is often described in terms of the Cosmological 
"Principle" , the application of either an FLRW or an LTB solution to the real Universe 
is better seen as a heuristic calculational strategy rather than a physical principle, with 
the risk of averaging-related artefacts occurring in both cases. LTB solutions provide 
an intermediate step between the FLRW solutions that force full homogeneity and an 
unknown family of generic solutions. 

Thus, here we consider LTB solutions. We first examine an LTB fit to recent 
observational data to see what it implies for the connectedness of spatial sections at early 
times (Sect. El). In Sect. HI we generalise this to a wider class of inhomogeneous solutions 
that contains two subclasses of distinct types of topological evolution (Theorem [I]) and 
present a conjecture and a corollary regarding one of the subclasses. We give case 

\ () denotes an implicitly denned average, () denotes functional dependence. 
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examples in Sect. 14.11 and Sect. I4.2[ using an LTB solution found previously [19], to 
show that the two subclasses of Theorem [1] are non-empty. The second subclass is 
considered in Sect. El Conclusions are given in Sect. El Unless otherwise stated, we only 
consider relativistic (non-quantum), comoving, dust solutions with a zero cosmological 
constant. 



2. Observational estimate 



There are many different fits of LTB models to observations — see Ref. [20] for a list of 
direct and inverse fits. Here, we consider a recent paper [l2| that phenomenologically 
used the inverse method to find an LTB model. That is, the authors started with 
functions implied by the FLRW model with Concordance Model values of the metric 
parameters [6[ and inferred LTB functions. A similar method and result are given in 



Ref. [21]. By construction, the two fits found in Ref. [12] provide good fits to the observed 
supernovae type la angular-diameter-distance-redshift relation and to an observational 
estimate of Hubble parameter evolution with redshift, based on differential stellar ages 
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of the oldest passively evolving galaxies at different redshifts 

The LTB models are comoving, dust models that, in general, allow the age of 
the universe in a given universe model in a comoving spatial section to be a function 
t — t B (r) that varies with the radial coordinate r. Thus, unsurprisingly, since the authors 



deliberately aimed to avoid making arbitrary assumptions, Figs 3 and 12 of Ref. [12 
show that the t^ir) solutions are not constant. In a comoving section at the present 
epoch to, the age of the universe increases from to at the observer to ~ to + 2 Gyr at an 
areal distance of about 3.7 Gpc. 

What are the topological properties of this solution? At times t > 0, let us assume 
that (i) the spatial section of the solution is simply connected, since E(r) > over the 



region of r > studied (Figs 2, 11 of [12J). Also, let us extend the solution by assuming 
that (ii) E(r) > Vr > 0. Thus, spatial sections at t > are the 3-manifold H 3 , with 
non-constant curvature. 



Figures 3 and 12 of Ref. 12] show that when —2 Gyr <, t < 0, a spatial section of 



the universe has a hole in the centre, where space has not yet emerged from the initial 
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singularity. For example, consider a spatial section at t — — 1 Gyr in Fig. 3 of Ref. 
In comoving coordinates, the closed 3-dimensional ball 

V = {(r,9,(f>) : r < r min (t = -1 Gyr)}, (1) 

where r m ; n (t = —1 Gyr) w 1.7 Gpc, consists of the initial singularity and a region of 
coordinate space beyond (earlier than) the singularity. The metric is only Lorentzian 
for t > tfi(r), i.e. r > r min (t), so the universe at t — — 1 Gyr is H 3 \V, i.e. a 3-manifold 
with a hole created by removing V from H 3 . 

Thus, this universe model evolves from H 3 \V to H 3 at early times. What is the 
areajjj radius R(t,r) = on the boundary <9V? This is given by (4), (7), and (2) [and 

§ The area of a zero thickness shell of radius r is 47r i? 2 (i, r) in (1) of Ref. thus, areal. 
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(6) for r = 0] of Ref. [12j|. As T](t,r) -»■ + , we have (j)(t,r) -> + and £(t,r) -> + , and 
thus -R(t,r) -)■ + , and t - t B (r) -» + , since E'(r) > and M(r) in Figs. 2, 4, 5, and 
6 are non-zero (for r > 0) functions of r only. Thus, the spatial volume of a shell at r 
shrinks to zero as t — > t B (r) + for fixed r, or as r — >■ r min (t) + at a fixed i. Within the 
spatial section H 3 \V at t, the boundary <9V appears as a single point. 

Relativistically, there is no problem with this solution. The high-redshift universe 
is born first, with the point r min (t) representing the unfinished early big bang process. 
The flexibility of the areal radius R(t, r) in LTB solutions allows comoving space to 
continuously be born from this singularity, which moves to successively lower values of 
r m ; n (t) as t increases up to t = 0. At t — 0, we have r m ; n (0) = and the singularity is 
replaced by ordinary spacetime points at (t > 0, r = 0). We can summarise this as 

r mi n(t) > 0, Vt < 
dr min (t)/dt < 0, Vt < 
V(t):=V[r<r min (t)], Vt < 

I d S 2 oc J R 2 (t,r min ) = 0, Vt<0, (2) 

JdV(t) 

where the pre-big-bang universe is considered to be non-physical and the metric is given 
in (1) of Ref. 12J. Thus, this universe model evolves from if 3 \{0} = S 2 x RtJ at 



t < to H 3 at t > 0. Comoving space is continuously born from the singularity until 
t = when the singularity disappears in the same way that it disappeared in parts 
of comoving space that were born earlier. Dropping the simplifying assumptions (i) 
and (ii) above does not make it possible to avoid the topology change in the former 
interpretation, since it just replaces if 3 by a more generic 3-manifold A4. 

3. Gaussian t-B{r,6,<f>) distribution 

The LTB solution presented in Ref. [12] (and |2l| ) is intended to demonstrate an 
example solution that fits key cosmological observations, but is not intended as a 
definitive replacement for the FLRW model with Concordance Model metric parameter 
values. Moreover, the LTB model is not a generic inhomogeneous model. An LTB 
solution constrained by more observational data could also be expected to have a more 
complicated non-constant £b function (unless this is imposed by assumption), while 
a more generic, inhomogeneous, dust solution could also be expected to have a non- 
constant £b function, as a function of three spatial variables rather than just one. While 



the solution [12J has only one (continuous) comoving spatial region where t-Q < maxiB; 
it is not obvious that this a requirement either of LTB models, or of more general 
cosmological (comoving dust) solutions of the Einstein equations. 

In solution [12j, lower densit y p tends to correlate with older regions of the universe, 



i.e. more negative £b (cf Fig. 3 of |l2j and the solid curves in Fig. 10 of [12]). A qualitative 



way to interpret this in terms of FLRW models is that for the same Hubble constant 
|| Defined here as R + := {x : < x <E R}. 



Inhomogeneity and topology 5 

H , a lower matter density fi m universe is older than a higher matter density universe 
This is only a qualitative guide to the LTB case, since both density and any typically 
defined equivalent of the Hubble parameter vary with t and r differently to the FLRW 
case. The same Figs 3 and 10 in 12] show that this qualitative inference does not always 
hold: lower p does not always correlate with more negative £b- 

To generalise from the solution [12| . let us suppose that a comoving dust solution 
to the Einstein equations has t-B(r,9,(p) drawn from a Gaussian distribution G(0,a), 
i.e. of mean zero and standard deviation a when smoothed on a length scale Ax. 
Gaussian density fluctuations on an FLRW background are a standard ingredient of 
modern cosmology, so even if it is unlikely that a given solution has a t-Q-p relation 
that is a function £b(p) (let alone a monotonic function), a Gaussian t-Q distribution is a 
heuristically reasonable hypothesis. Now consider an approximately flat, cubical, small 
region of side length 3 Ax of which the central (Ax) 3 small cube contains a region with 
t-Q < —3a, i.e. born unusually early. The probability that this small cube is connected 
to another small cube with ts < —3a, i.e. that it is not isolated by iso-bang (constant 
£b) contours, is the complement of the probability that the 26 small cubes around it all 
have t B > -3a, i.e. P = 1 - Q [l + erf(3/ v / 2)]) 26 « 3%Q The chances that the second 
t B < —3a cube touches a third cube outside of the original (3 Ax) 3 region, and that the 
(n > 2)-th cube touches another small cube yet further away for n > 3, rapidly decrease 
with increasing n. 

Thus, te = — 3cr iso-bang contours will tend to form isolated 2-surfaces. That 
is, with a fixed smoothing scale Ax and in a large enough comoving spatial region, a 
Gaussian distribution in t-Q implies that there will tend to (statistically) exist many 
regions (3-volumes) {V{\ with t B < — 3a that are spatially isolated from one another 
within the 3- manifold. At t ^> 0, we label the latter Ai. Let us assume that Ai is 
connected and that its volume is ^> (Ax) 3 . 

Now consider the comoving spatial section at t = —3a. The boundaries of the 
regions {V{\ with t-Q < —3a {dVi} are 2-spatial iso-bang contours. These regions have 
already emerged from the initial singularity since t—t-Q > 0, apart from their boundaries, 
where t — t-Q = 0, i.e. at the initial singularity. Since the {Vj} are isolated from one 
another, they constitute a set of disconnected 3-manifolds. Hence, the universe at 
t = —3a consists of the spatially disconnected 3-manifold U{Vj}. 

The choice of — 3a is for illustration only. Any reasonably high x ,> 3 will 
(statistically) give a spatially disconnected universe at t — —xa, given a large enough 
spatial volume and a Gaussian distribution of £b as stated above. At the same time t, the 
parts of space U {Vj} have not yet emerged from the initial singularity. If we follow 
the spatial section back in time from t = —3a, then the boundaries dVi correspond 
to t = —xa for increasing x, i.e. they shrink smoothly, possibly subdividing further, 
eventually vanishing into the singularity. For t < minted, 0, </>), the global bang time, 
no more V, exist. 



For a fixed cosmological constant Q\. 
+ Conservatively, cubes that touch at corners are considered to be connected. 
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Moving forward in time, how do the V, merge together? The boundary dVi for 
any disconnected region has zero 2-surface area, as in the case of R(t, r min ) = in the 
solution 12]. That is, the boundary of Vj is S 2 with zero 2-surface area, i.e. a point-like 
singularity. Thus, Vj can be thought of as a 3-manifold with one point excluded. For 
intuitive purposes, it can be useful to think of an azimuthal equidistant projection of 
the Earth's surface, centred at an arbitrary geographical location, with the antipode 
corresponding to the big bang initial singularity. The antipode can be thought of either 
as a large, coordinate-space, zero- circumference circle that bounds the 2-manifold from 
the "outside" in the projected map, or as a single missing point "on" our usual intuition 
of the Earth's surface. 

Again, as in the solution 12|, comoving space is born from this singularity, so that 
there is comoving growth of the spatial region Vj. As t increases, the threshold for 
the iso-bang contours increases, so that the Vj eventually touch and pairs of Vj merge 
together. Again writing t = —xa, as — x becomes more positive, t reaches a high enough 
— xo such that the probability for an isolated = —xa region to exist becomes 
negligible and the universe becomes fully connected. 

For a Gaussian £b distribution, it could be expected that the zero 2-surface area 
of the boundary of an isolated region V» at t = —xa will tend to be topologically S 2 
in coordinate space, so that Vj is topologically S 3 \dVi = S 3 \{0} However, other 2- 
manifolds for the coordinate space representation of dVi are possible, e.g. a 2-torus, so 
that Vi is S l x (S 2 \{0}). 

Thus, we find that if a universe described by the Einstein equations is £b- 
inhomogeneous, then, even with several simplifying assumptions (Gaussian distribution 
of t-Q at a given smoothing length, Ai connected and simply connected for t ^> 0), 
there is a very high probability that it emerged from the (spacetime-smooth) mergers 
of comoving spatial sections that were spatially disconnected from each other prior 
to their mergers. This is distinct from questions of causal connectivity and raises 
fundamental questions for the physical properties of the spacetime model that are 
considered acceptable. 

Let us propose a generalisation. 



4. Relativistic, post— quantum-epoch topological evolution 

Theorem 1 Let us define the generic class A where {g~} G A if is a (regular) 
extension to t-g < t < to of a dust (pressureless) metric g\t , i.e. g~ solves the Einstein 
equations over t B < t < t , where g\ to on a comoving spatial section (3-manifold) at 
~ t , which we call M., 

1. solves the Einstein equations, 

2. is regular, and 

3. has an approximately homogeneous density p. 
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Here, t is the age of the Universe at the location of our Galaxy, and t B is a function 
of comoving spatial position defined by the initial big bang singularity. Two non-empty, 
distinct subclasses of A are A d and A m as follows. 

(i) A d , in which the universe is born from an initial singularity at t — )■ (mint B ) + 
as two or more spatially disconnected regions (3-manifolds) Vi, each of which is 
bounded by at least one singularity ( of zero spatial volume ) from which comoving 
space emerges continuously. The Vi successively merge together to form the t > 
max£ B > mint B connected 3-manifold Ai . The Vi themselves are born, in general, 
at different times, and their enumeration changes as a function of t, because of 
their mergers. 

(ii) A m , in which the universe is born from an initial singularity at t — >■ (mint B ) + 
as a connected, simply connected region V\ bounded by at least two singularities 
during mints < t < mint B + St for some 5t > 0. Comoving space is continuously 
born from the singularities, which join together smoothly in pairs (or n-tuples, with 
n > 2), so that the spatial section at t > maxte is a connected, multiply connected 
3-manifold M. . 



The heuristic Gaussian t B discussion suggests that (i) of Theorem [T] is correct, 
without establishing it rigorously in an exact solution of the Einstein equations. Neither 
(i) nor (ii) of Theorem [1] are relativistically problematic. However, both (i) and (ii), 
if they are correct, imply a new type of causality problem distinct from the particle 
horizon problem: teleology. How is it possible for the growing singularities in spatially 
disconnected regions to "know" where other regions and their singularities are "located" 
in order for the singularities to join together? A coordinate system such as that used 
for LTB models is convenient to work with, but if there are spatial islands in the spatial 
part of the coordinate system, the remaining "sea" of space consists of a purely fictional 
construct — useful for coordinate-based intuition — until comoving space is born there, 
converting it from Active, coordinate space to physical (metric) space. If we only have a 
relativistic spacetime (with a Lorentzian metric everywhere), then individual Vi cannot 
"know" that they must be embedded in a future coordinate system that will unite them. 
The transition from simple to multiple connectedness (ii) is physically simpler, since the 
singularities exist in the same, connected initial manifold. However, a physical theory 
is still required in order for the singularities to "know" where their future partners are 
located in the same space. Suppose a given function t B varies by within a "few" orders 
of magnitude of the variation in the empirically derived t B m the solution 12] . If a "few" 
means <C 10 60 , then the theory would have to allow singularities across spatial regions 
that are significantly post-Planck to connect to one another. For a "few" meaning 
<C 10 10 , the epoch would be significantly post-nucleosynthesis. 

Examples of members of A d and A m are given in Sect. 14.11 and Sect. 14. 2\ proving 
Theorem [TJ This provides the basis for hypothesising that solutions with non-constant 
t-Q are more common than those with constant t B . 
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Conjecture 1 For a measure fx on A that is physically motivated at late epochs (and 
that does not contradict early disconnectedness) , the measure of solutions that are not 
primordially disconnected is small, i.e. fi(A\A A ) <C fi>(A). 

Corollary 1 If 

(i) Conjecture Q is correct, and 

(ii) g for our real Universe is chosen randomly from A, and 

(Hi) the variance of the t-Q timescale is 3> 1CT 60 times that estimated empirically in 



Ref. \l£], 



then spatial disconnectedness occurred at early epochs, and this happened at an epoch 
much later than the quantum epoch. Thus, the inverse method of using extragalactic, 
astronomical observations to extrapolate back towards the initial singularity, for example, 
numerically using the (3+1) -formalism [e.g. 25], would be likely to yield evidence of 
early, post-quantum spatial disconnectedness . 

However, as commented above, mergers of disconnected "universes" require either 
(i) teleology or (ii) a physical theory beyond general relativity. Hence, if we reject 
(i) teleology as unphysical, and if the conditions of Corollary CD are satisfied, then 
the relativistic extrapolation of astronomical observations towards the past would 
be expected to provide observational evidence of a spacetime theory beyond general 
relativity [not necessarily an (n > 4)-manifold] at a post-quantum epoch. That is, 
inhomogeneous models inferred from observations would typically find pairs of regions 
of the observed universe (e.g. the cosmic microwave background) to be spatially 
disconnected at early epochs, not just causally separated. 

It also follows that if teleology is rejected by assuming constant £b, but the three 
conditions of Corollary [1] are satisfied, then the constant t^ assumption would force 
a wrong solution to emerge from extrapolation of observations into the past. It would 
prevent the correct inference of a post-quantum-epoch spacetime theory beyond general 
relativity from being made. 

Conversely, observational limits on the existence of (primordial) singularities and 
inhomogeneities surrounding them [e.g. in the observed cosmic microwave background 
(CMB)] should be used to constrain the allowed variance in £b for observationally 
realistic inhomogeneous cosmological solutions of the Einstein equations. Since the 



solution [12| has maxt B ~ mint B ~2 Gyr over about 4 Gpc, several singularities 
corresponding to the initial big bang singularity should exist in the observed CMB, 
which has a comoving surface area of about 47r(10/i _1 Gpc) 2 and an age of about 
0.0003 h^ 1 Gyr according to the FLRW phenomenological fit, if this particular solution 
were correct. For consistency, the calculations would have to be performed with a 
background-free approach (i.e. without assuming an FLRW model). 
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Figure 1. Example of iB-inhomogeneous positively curved LTB solution (1^1 ITU]) that 
is born as disconnected spatial sections that smoothly merge together via the initial 
singularity (cf Fig. 7(a) of [l^]), showing a finite part of the comoving spatial section, 
which is of infinite length in the r direction. The universe exists (has emerged from the 
big bang singularity) in the shaded region (excluding the singularity itself, appearing 
as a sinusoid here). The thin horizontal line shows a spatial section of the universe at 
t = — 0.42/i _1 Gyr, during which some parts of the universe exist, and others do not 
yet exist. 
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Figure 2. Density p (left) and a Hubble-like parameter \H\ (TH21) (right) of the 
solution shown in Fig. [T] Values are plotted for the ranges of r where the universe 
exists, i.e. t > £b(?")- Values increasing arbitrarily are shown to limits that depend on 
numerical implementation details and avoid obscuring the legends. A late epoch, near 
recollapse, is shown by a thinner curve in both panels in this and the following figure; 
its H parameter is negative. 



4-1. Spatially disconnected sections that merge 

We show that class A d , i.e. (i) in Theorem [1] is non-empty, using an explicit example 



of the "string of beads" LTB solution [19| (see also [24j). This is a positively curved 



solution. This class of solution requires the radial metric component g rr to be defined 
as a limit, because of behaviour at what (in the FLRW case) is the model's equator [e.g. 



251 ] . This particular example has a function with sinusoidal behaviour, with all the 
minima and maxima occurring at a single pair of values, mints and maxte, respectively. 
This is not a general requirement, it is just a characteristic of this particularly simple 
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solution. 

Let the metric be 



ds u 



-dt 2 + lim 



R 2 r (t,r) 



-dr 2 + R 2 (t, r) {d9 2 + cos 2 6 d(j> 2 



(3) 



>r 1 + 2E(r 

1, the gravitational constant G is written explicitly and E(r) is a curvature- 



where c 

related function, e.g. (1), (21), (23) [17J; (2.1) (26|; (1) (12 
Einstein equations exists if 

2GM 



1) [25|. A solution to the 



p 



2E + 

M, 



i? 



4tT.R 2 .R , 



where 



i2(*,r) 
*s(r) 



GM(r) 
' 2E(r) 
GM(r) 



[1 — cos7](t, r)] 

■e(«,r) 



(4) 
(5) 

(6) 



(7) 
(8) 



[-2E(r)] 3 / 2 
r)(t,r) — smrj(t,r 

where M(r') is a weighted integral of the density over < r < r', tB( r ) is the bang time 
and ^ and r] are auxiliary functions. 
Following Section 8 of [l9j], define 

1 



E(r) 
M(r) 

M n 



■[1 -Ei sin 2 (r)] 
M (l + Mxcosr) 
F ^ + GM (1-M 1 ; 



2Gi/o(fi m - 1) L5 (9) 
where the FLRW dimensionless matter density parameter Q m := 8tiGpo/(3Hq) is set 
to l] m = 2, po is the present matter density, Ho is the FLRW Hubble constant, and 
Mq is chosen to get a time scale roughly comparable to that of an FLRW positively 
curved model with zero cosmo logical constant. In order to avoid R r having zeroes 
where M r does not have zeroes (see Section XIV. B of 27]), the second derivative of 
(l + Micosr)/(l-£isin 2 r) at must be negative, i.e. E x < 0.5Afi/(l +M X ). Thus, to 
obtain a sub-Gyr time scale of variation in £b ; i-e. comparable to (but more conservative 
than) the solution 12], the parameters are set at 

Mi = 5 x 10~ 5 

E l = 3 x 10~ 6 . (10) 

Figures [U [2] and [3] show this solution at early epochs and the evolution of some key 
properties. The 3-Ricci scalar is 



RR , 



+ 



-) 

R 2 J 
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Figure 3. As for Fig. [2 the Ricci scalar 3 i?/6 (fTTj) (top-left), the areal radius 
R(t,r) (top-right), and the radial metric component y/g rr (t, r) (bottom- left). The 
radial proper length d(t, r) (Tl3|) at several pre- and post-connection epochs t is also 
shown (bottom-right; the thick curves match two of the epochs shown in the other 
panels, as labelled). 



and following 271 ]. a Hubble-like expansion parameter is defined (14), (29) 27 



3 \ R R r J 

The early epoch curve in the top-right panel of Fig. [3], i.e. for R{— 0.42/i _1 Gyr, r) 
shows numerically what can be seen in (jGj) , (GO), and (El): provided that the factors that 
include E and M are well-behaved, the one-sided limit rj — > + <^ £ — > + corresponds 
to R(t,r) — > + , and t — t B (r) — > + . Thus, as for the solution [12J, zero-surface area 
2-spheres, i.e. point-like singularities, bound the post-big-bang parts of the universe 
model. 

In coordinate space, it is clear that the spatial sections of the universe are 
disconnected at t < 0. What happens at and near the coordinate points (t, r) = 
(0, (2n+ l)vr), n G Z? Let us, w.l.o.g., consider (t, r) = (0, 7r). Since R = at this point, 
the metric ([3]) has a non-Lorentzian signature: this is the initial big bang singularity 
from which the comoving spatial point (t > 0, tt) is born. Thus, at t — 0, the two parts 
of the spatial section (0, — it < r < it) and (0, it < r < 2tt) are disconnected from one 
another by the point (0, 7r) in coordinate space. While disconnection by a single missing 
point might seem trivial, since mathematically, adding a point to a manifold can remove 
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Figure 4. As for Fig. [TJ example of inhomogeneous positively curved LTB solution 
with non-simultaneous (i.e. non-constant) tB{r) that is born with a simply connected 
spatial section that smoothly connects to itself via the initial singularity, becoming 
multiply connected (S 2 x 5 1 ), shown in the universal covering space. One copy of the 
fundamental domain (FD) lies between the two vertical lines. 



a singularity, the physical significance would be non-trivial. The addition of a single 
point "at infinity" to infinite Euclidean 3-space M 3 is enough to transform the latter 
into S 3 , although physically, this would be absurd. 

How does the radial component of the metric behave near the connection points 
(0, (2n + l)vr)? For t > 0, Fig. [3] shows the anisotropic way in which the metric evolves. 
As t — > + , R decreases (top-right) but g rr increases (bottom- left). The latter increases 
without bound as t — > + at (t, (2n + 1)tt) and remains infinite at the r boundaries of 
the disconnected sections. However, the integrated proper length 



from the centre of an initially disconnected section to its boundary, i.e. to the big-bang 
singularity, remains finite (bottom-right). 

Unless a pre-big-bang scenario is introduced, the comoving spatial sections of the 
universe during mint B < t < consist of the disjoint union U ig z S 2 x (0,1). This is 
not an issue of particle horizons within acausal spatial sections; the spatial sections are 
disconnected. With the parameters chosen, the delay before these grow and merge with 
the "rest" of the future-to-be-created spatial section is more than 100 h~ x Myr, i.e. long 
after nucleosynthesis. 

Thus, A d of Theorem [T] is a non-empty set. The age of the Universe to used in this 
example is that for an FLRW model with Q m = 1.015, = 0; a to =10 h~ l Gyr model 
can be calculated trivially by modifying M . 




(13) 



4-2. A connected, simply connected section that becomes multiply connected 
Taking the solution ([9]), (JTDjl . apply the holonomy 
y:(t,r,9,<f>)^(t,r + 2ir,e,<f>). 



(14) 
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The spatial sections for t > are multiply connected, i.e. S 2 x S 1 . This is an exact, 
non- vacuum solution of the Einstein equations with a multi ply connected spatial section, 
similar for t > to the S 2 x S 1 solution published earlier |28|l*l 

But at t < 0, the spatial section is S 2 x (0, 1) i.e. it is a single, connected, simply 
connected 3-manifold. Hence, a simply connected universe can smoothly become multiply 
connected at early (post-quantum) epochs: the class of solutions A m of Theorem [1] is non- 
empty. Figure 0] shows the universal covering space of this solution. 

5. Evolution of a connected 3-manifold 

Let us consider a more conservative hypothesis than Conjectured], i.e. a hypothesis that 
rejects primordial disconnectedness as unlikely, but does not force t# to be constant. 

Conjecture 2 For a measure /i on A that is physically motivated (and does not 
contradict the mergers of early epoch singularities ), 

(i) the measure of solutions that are disconnected is zero, i.e. fi{A d ) = 0, and 

(ii) the measure of the class of solutions with constant ts on comoving (always 
connected) spatial sections is small, i.e. /i(A\A tB ) <C fi(A tB ), where A tB is the class 
of solutions with non- constant t-& and spatial sections that are always connected over 
min £b < t < t . 

By the definition of A m (Theorem HJ, A m C A tB . 

If Conjecture [2] is correct, then a universe is most likely to be born connected at 
mints, with at least one singularity that disappears later (as in Sect. [2]). If the universe 
is born with many singularities, then some may disappear individually (as in Sect. |2]), 
some may disappear in pairs (as in Sect. I4.2jl . and others could, in principle, disappear 
in n-tuples with n > 2, even though it may be hard to find exact metrics as examples of 
regular mergers of n > 2 primordial singularities. Thus, if Conjecture |2] is correct and 
if the universe is born with many (N ^> 1) singularities, then an example of a minimal 
extension of general relativity that would describe the evolution of the universe would 
be a definition Vi, j G Z : i, j < N of 

(i) Pl(g(t), t), the probability that singularity i at time i disappears at time t in a way 
such that g is regular W < t + 5 for some 5 > over the whole spatial section, and 

(ii) \fn : 2 < n < N, P l n 1 ' i2 '-' l "(g(t), t), the probability that the singularities i\, i 2 , ■ ■ ■ , i n 
at time i merge together at time t in a way such that g is regular W < t + 5 for 
some 5 > over the whole spatial section. 

Given the existence of the numerical solution in Sect. [2]) and the analytical solution in 
Sect. 14.21 and the requirement that in the latter case, the two pre-merger metrics must be 
post-merger compatible, it would seem reasonable that P\ ^> Pi and % < j =^ Pi^> Pj, 

* The solution in Sect. VI. D of [28j], with E non-differentiable at d/A and 3d/4, is stated by the 
authors to be a "3-torus", presumably by analogy with S 1 x S 1 =: T 2 . However, the analogy fails, 
since S 2 x S 1 ^ T 3 := ^x^x S 1 . 



Inhomogeneity and topology 



14 



although this is only speculation. Two obvious classes of models would be those that 
define the probabilities Pi and P n ,n > 2 to be independent of the (comoving) spatial 
locations of neighbourhoods of the singularities, and those that define the probablities 
to be dependent on the spatial locations or on global properties (e.g. mean 3-Ricci 
scalar, topology) of the spatial section. The P\ and P n , n > 2 could also depend on the 
comoving spatial number density of the singularities. 

A model of the functions P n ,n > 1 would provide a minimal extension of general 
relativity which could be used to calculate the probabilities that a universe evolved 
from a connected, simply connected spatial section to a connected, multiply connected 
spatial section, and which topologies would be mostly likely to remain at t > maxte- 
If the P2 and the number of spatial singularities are high enough, then evolution 
to a multiply connected spatial section would become more likely than to a simply 
connected spatial section. If the variance of the t-Q timescale is ^> 10 -60 times that 



estimated empirically in Ref. [12j, then this model would apply at significantly post- 
quantum epochs. Nevertheless, the requirement of discreteness (since the singularities 
are discrete within comoving spatial sections) and the suggested probabilistic nature of 
the model suggest a quantum model, for example, with a Heisenberg-like singularity 
merger principle 

AH(i,j)As(iJ)>h, (15) 

where AH(i,j) is the absolute difference between the Hubble-like parameter (fl2l) 
averaged over some comoving neighbourhoods of singularities i and j, and As(i,j) is 
the shortest spatial, comoving geodesic distance^! between singularities i and j. 

6. Conclusion 

We have examined the early epoch topological evolution that corresponds to non- 
simultaneous big bang times in non-empty, inhomogeneous dust models of the Universe 
using a recent empirical estimate and an older analytical exact solution. Even if 
the variance of t-Q estimated empirically is overestimated by several tens of orders of 
magnitude (in their introduction, the authors suggest that a more realistic time scale 



would be ~ 100 yr, i.e. about 10 7 times smaller than their empirical solution 12] ), it is 
still post-quantum unless £b is constant to within the Planck time scale, i.e. ~ 10~ 43 s. 



Other work also estimates a Gyr time scale [e.g. Figs 6, 8, 29]. Thus, the topological 
evolution implied by te-inhomogeneous models has fundamental physical implications, 
which in one case (teleology, for evolution from a disconnected to a connected spatial 
section without extending the spacetime model beyond 4-dimensional general relativity) 
violates Lorentzian causality. 

The physical implications could be avoided by imposing £b = constant as an 
assumption in cosmological modelling. One problem in assuming constant t-Q is that for 
flat LTB solutions, generalisations beyond the FLRW model are rejected. That is, the 

jj Two points in a 3-manifold may be joined by more than one distinct geodesic. 
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combination of E(r) = and t^ir) = leaves no freedom to adjust the third "arbitrary" 



function M(r); see VIII (63a), XIV.B in Ref. j27j. Section XIV.B of Ref. (27| also 
discusses the restrictions on LTB models implied by imposing £b(?") = Vr in the positive 
and negative E(r) cases. More importantly from a physical point of view, allowing 
te spatial dependence to be a result of comparison between models and observations 
rather than an assumption could potentially lead to evidence that extensions to general 
relativity are required, which would otherwise be artificially suppressed. 

We have formalised some of the possible properties of subclasses of solutions of 
this type and of possible implications in Theorem [1] Conjecture [IJ Corollary [1] and 
Conjecture El Conjecture [2] opens the way to calculations of the probabilities of a simply 
connected initial spatial section smoothly evolving into a multiply connected spatial 
section, based on a choice of functions P 1 4 (g(t), t), \/n : 2 < n < N, P^ 1,42, - '* n (g(t), t) 



as defined in Sect. E 



30|, |3l|, |32|, |33j, |3J, |35 



39 



above. This would have considerable observational interest [e.g. 

since it would offer an alternative to the topological 
for theoretical understanding of the topology of the 



36 



31 



40 



38 



41 



acceleration effect 
Universe. 

How was it possible that post-quantum-epoch topology change without causality 
problems was overlooked in cosmic topology literature? It has generally been thought 
that if the spatial sections of the Universe are compact, then the topology of spatial 
sections of the Universe cannot have evolved at post-quantum epochs, because this 
would imply the existence of closed timelike curves or a discontinuity in the choice 
of the forward light cone, as a consequence of Geroch's Theorem 2 42], [Sect. 9.4.1. 



43J, 44j, and both are generally considered unphysical. Singularities make spacelike 



sections non-compact, so that the theorem no longer applies, but it is not immediately 
obvious that an astrophysically realistic black (let alone white) hole could change the 
global topology of spatial sections in a way that leads to approximate homogeneity 
in the late-time Universe. What was overlooked was the fact that a non-constant 
£b provides (in general) an arbitrary number of singularities in early, post-quantum 
comoving spatial sections, which in spacetime constitute just one singularity — the initial 
big bang singularity which is generally accepted as physical in relativistic, non-quantum 
cosmology. Moreover, as illustrated above, the density and curvature inhomogeneities 
near vanished singularities/connection points can become much weaker, i.e. enter a 
"decaying mode" in a perturbed FLRW context. The LTB models provide a useful tool 
for studying examples of characteristics that are counterintuitive for FLRW-like models. 

Although in this work we only consider topology change implied by non-constant 
ts in classical relativity, see Ref. 44J for a quantum gravity approach using sums of 
histories and Morse theory. 
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